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IMPULSIVE MOTION ELASTO-PLASTIC 
H.H. Bleich M.G. Salvadori 


SYNOPSIS 


general analysis the vibrations elasto-plastic 
structures given and applied free-free beams moving under 
dynamic forces. The motion during the successive elastic and 
elasto-plastic stages expressed terms the normal modes 
the structure and the maximum plastic angle evaluated, 
The case impulsive loads (initial velocities), which 
practical interest, specifically considered, and means 
simplifying assumptions, upper and, often, lower bounds the 


Plastic are 
GENERAL THEORY. 


Consider ideally elasto-plastic body capable vibrat- 
ing elastically single infinity normalized, orthogonal, 
subjected system external time varying forces, symboli- 


cally indicated F(t), the vectorial elastic displacement 


one its points may expanded into series 


where the generalized coordinates the absence damping, 


modes 


satisfy the Lagrangian equations 


This paper presents the results research conducted 
Columbia University under contract Nonr 266(08) with the 
Office Naval Research. 


Professors Civil Engineering, Columbia University, N.Y. 
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Eq. (2), the kinetic energy and the potential 
energy the elastic body, and the are the generalized 
forces corresponding the system external forces 

Let now assume that time plastic state 
the point while all other points the stresses remain 
elastic. Then for may consider the constant force 
which the continuity destroyed the area far the 
normalized modes this new, relaxed body its 
generalized coordinates and the generalized force 


corresponding F(t) P,, the vectorial elasto-plastic 


4? 
displacement may expanded series 

oo 

(3) 


where the satisfy the Lagrangian equations 


The continuity the displacement and the velocity 
require that: 


The modes the relaxed body may conveniently 


obtained, many cases, the method "complementary 
systems" (see H.H. Bleich, "Frequency Analysis Beam 
and Floors", Transactions, ASCE, Vol. 115, 1950, 
1023 
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Multiplying Eqs. (6) integrating over the entire mass 
the body and remembering the orthogonality conditions the 
modes 


the initial conditions for the reduce to: 


where: 


Similar relaxations will give the equations motion when 
other areas the body become successively plastic, provided 
the plastic areas not change shape location, 

The elasto-plastic motion defined Eq. (3) will continue 
until the stresses the plastic area back into the 
elastic range. this takes place time, say 
the displacement for >t, 
new generalized coordinates satisfying Eqs. (2) and the 
initial conditions continuity 


1 
given Eq. (1) with 


Co Co 


Multiplying these equations and integrating over the 
whole mass the body, the initial conditions for the 
become, because the orthogonality the modes 


(8) 
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i= 
oo 


(12) 


sequence changes from the elastic the plastic range and 
vice versa allows the description the elasto-plastic forced 
vibrations the body any time 

Under particular circumstances, one more the plastic 
regions may either change shape location both before 
reversal stress before additional plastic regions appear. 
these cases (which are not considered the present paper), 
new and more complicated equations motion must set and 
their integration may have performed approximate, step- 


by-step methods. 


THE ELASTO-PLASTIC VIBRATIONS FREE-FREE BEAMS. 


typical example application the general theory 
Sec. consider the forced vibrations free-free beam 
length variable flexural rigidity and variable mass 
per unit length (Fig. 2). 

The free vibrations such beam are governed the 
differential system: 


where primes denote differentiation with respect the 
variable measured along the centroidal axis the beam with 
origin its middle, and dots denote differentiation with 


respect time. 


The substitution 


(11) 
where 
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Eqs. (13), (14), and (15) shows that the beam can vibrate 


infinite number modes satisfying the equations 


and the boundary 


The modes are orthogonal and may normalized that 


convenient indicate the subscript the 
number nodes the i-th mode; thus, represents the 
rigid translation the beam, its rigid rotation, 
its first elastic mode, etc. 

Indicating the general integral Eq. (17) 

and substituting this value Eqs. (18), obtain 
homogeneous system linear algebraic equations the con- 
stants whose determinant set equal zero determin- 
the frequency the i-th 

The modes can evaluated rigorously from the diff- 
erential system (17), (18), this impossible in- 
convenient, they can obtained approximate 
experimental methods. 

beam due time varying pressure p(x, per unit 


modes 


Substitution (21) the Ligrangien (2) and the 
initial (15) shows that the gem ralized coordinates 


satisfy the 
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given the integral: 


Multiplying Eqs. (23) and integrating over the beam 
the conditions for the become, Eqs. (19): 


(25) 


any time during the elastic stage the vibrations, 

Under the assumption ideal elasto-plastic material 
the elastic stage the vibrations ends the smallest 
the time say which the moment reaches the 
plastic capacity value anywhere the beam. the moment 
reaches the capacity value first section the 


clastic stage the vibrations ends when 


the capacity moment the beam the section 
If, least for time, the capacity moment not reach- 


any other scction, the moment M(x, will remain constant 
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and equal The elasto-plastic vibrations the beam 
are then identical with the vibrations another beam equal 
the beam considered, but with hinge and acted upon 
the same load p(x, and two external, equal and opposite 
moments the hinge. 

Indicating the orthonormal modes the hinged 
beam, the vibrations the original beam may 


expanded into series 
oo 


where the new generalized coordinates satisfy the differ- 


ential 


being the frequency the i-th mode the hinged beam, 


and the generalized forces due the load and the moments 


The generalized force thus given 


The continuity displacements and 


Substitution Eqs. (21) and (28) these conditions, 
plication and integration over the whole beam, 


reduce these the 


t 
(42) 
wnere the are given Eqs. (9). 
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When beam has mass symmetrically distributed about its 
middle point and the pressure p(x, also symmetrical and 
such that the first hinge will develop the middle point, the 
modes the hinged beam are obtained noticing that 
there force transmitted through the hinge and that each 
portion the beam vibrates free-free beam length 
The modes are the modes the two half beams normalized 
over the 

While the moment remains constant, plastic angle 
develops this section; its value equals the dis- 
continuity the slope the elasto-plastic deflection 


(33) 


When time, say this angle reaches maximum value, the 
relative motion the two portions the beam reversed and 
the moment becomes elastic Under the assump- 
tion small plastic deflections and elastic stresses 
all sections the beam except the beam deflections 
are given again the series (21) with new generalized coordi- 
nates satisfying Eqs. (22) and conditions continuity 
cannot carried out with more than the first few terms the 
series expansions (21) and (28) terms are usually 
sufficient), conditions continuity the displacements 
and the velocities not necessarily make the 
accelerations continuous that time and hence may give 
numerical discrepancy the values the moment M(x, 
computed inertia forces the elasto-plastic range and 
according Eq. (26) the elastic range. therefore 


often advisable obtain the initial values the 
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imposing conditions continuity the and the 


accelerations, order obtain continuous 


ard, means (21), (28) and (19) reduce to: 
os 


120 


Under the assumptions the stresses the beam may 
evaluated any time terms successive elastic and 
elasto-plastic the type given Eqs. (21) and 
(28). some cases the successive development plastic 
hinges more then one section must considered, but the new 
hinges can treated shown above and introduce essential 
difficulty the solution. cther cases, one more the 
successive hinges may shift their locaticn continuously that 
whole portion the beam becomes these 
ing hinges difficulties the integration 
the which then can e.g. solved step-by-step, 
using integration Alternatively, such cases 
can assuming that beam has one 
more points the plastic region; all plestic de- 
formations then points, and moments 


slighly excess the points are tolerated. 
(1) Cee, Lee and Symonds. Plastic 
Deformations wnder Transverse 
Salvadori and Diliaggio, "On the 
Plastic Hinges Beams Quart. 
App. 1953, 223 ff. 


(35) 
where 
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MOTION FREE-FREE BEAMS WITH CONSTANT FLEXURAL 
The stage. 
impulsive motion beam due forces high in- 
short duration equivalent, the limit, the 
the previous sections shall derive 
governing the elasto-plastic motion free- 


given the symmetrical distribution 


modes such beam are given 


rt 


the roots the transcendental equation 


Cc., Inc., New York, second edition, 


(Fig. 
cos cosh k,x 
sin k,x sinh k,x 
and the frequencies by: 
TABLE 
| 
ig 
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cos cosh (40) 
The coordinates the absence the 


external p(x, satisfy the homogeneous equations 


and, since assumed equal zero, the initial condit- 


2r) 


Hence, the are given by: 


The first elasto-plastic stage. 

The symmetrical initial velocity (37) excites only the 
even modes the elastic beam and creates bending moment 
neighborhood Hence the section reaching plasticity 
first the section and the plasticity time" 

which the first plastic hinge occurs determined the 


equation 


The modes the hinged beam are noticing 
that each half-beam vibrates like free-free beam and that 
their common end must have the same deflection. Hence each 
half-beam has modes obtainable from those Eqs. (38) 


The roots Eq. (40) approach asymptotically the values 


fe 
\ 
{ 
(0) 
(42) 
(0) 
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changing into Thus, the referred the 


middle point nalf-beams, 


(45) 


The elasto-plastic generalized coordinates satisfy the 


non-homogeneous equations 


and the initial conditions (32), which the constants are 
computed substituting the values and Eqs. 
(38) and (45) Eq. (9). The generalized forces due the 


moments are given 


where the angle equals the slope discontinuity the i-th 


0 
ae 2 ' 
with 
= 
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Hence the generalized coordinates take the form 


The elasto-plastic displacements the beam are computed 
substitution the given Eqs. (45), and given 
Eqs. (50), (51) Eq. (28) and the plastic angle 


computed from Eq. (33), which becomes this case 
>) 


The maximum the plastic angle 

practical importance both because determines the 
amount energy, Spent plastic flow and because 
the time, say its occurrence establishes the transition 
from the elasto-plastic the elastic The time 
determined setting equal zero the derivative 
given Eq. (52), that is, the equation: 

upper bound for may easily obtained assuming 


that the available energy impact which not transformed 


(50) 
where 
| 
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into energy goes entirely plastic flow. 
allowing that the energy would actually spent 
elastic vibrations for plastic work, the 


tained larger than the true The 


computed from Eq. (52) con- 


rigid rotation mode The plastic 
irs which time the total energy the 
l 


the beam having plastic hinge the center but otherwise 
rigid, che amount energy available for deformation will 

energy energy the second mode mist 


wh 


therefore give bound for the angle Thus with 


2 


The maximum absolute value the upper bound 
time: 
and equals: 


(57) 


max 


where given Eq. found from Eqs. (48) 
= 2 ' 


The good approximation the true 
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q 
upper 
(0 
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that the maximum bending moment which would develop the 
beam never became plastic times larger than the 
capacity moment. Thus even the impact time very short, 
the maximum plastic angle may well approximated the 
plastic approach, although the time history the in- 
correct due the assumption plasticity time 

Another approximation which gives more correct 
time history and often lower bound, may obtained 
the beam into rigid plastic body rather 
than this case obtained again from Eq. (52) 
considering only the rigid rotation mode from 


on: 


The absolute value becomes maximum time: 


and 


where and are given Eqs. (32). 


ing elastic moments several times larger than the capacity 
moment 


The successive elastic and elasto-plastic stages. 


From on, the deflections become elastic again and 
are computed means the elastic modes (38) and new 


generalized coordinates satisfying the initial conditions 
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(35), 
where 


(=A 
the are given Eqs. (51) and the (9). 


most cases impulsive the consideration 
further elasto-plastic stages unnecessary, because the 
energy spent plastic flow, but successive elasto-plastic and 
elastic stages may determined whenever needed equations 


Eqs. (32) and (35). 
EXAMPLE. 


The theory presented Sec. applied specific 
beam order illustrate its main features. The results 
this application, which the first terms the series were 
taken into account, are given graphically this section. 

The beam considered has the fcllowing characteristics 
consistent system units: 

Length: 320 

Mass per unit length: 
Capacity moment: 6.69 10" 
Fundamental period the beam: 1.36 


Two values the initial velocity parameter were used: 
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evaluated from Eq. (44), 


The initial plasticity time 


was found 


Fig. gives the graph the elastic moment fora 
full period the fundamental mode the beam, under the 
assumption infinitely large capacity moment throughout 
the Fig. gives the graph the elasto-plastic moment 
for 11.1 and for full period the fundamental 
elastic mode the beam. seen from Figs. and that 
the maximum elastic moment would about twice the cap- 
acity moment. 

The moment becomes elastic again time 0.375 given 
Eq. (53) and does not revert its plastic value within the 
time considered, 

Fig. gives the plastic angle given Eq. (52), 
together with the upper bound Eq. (55) and the bound 
Eq. (58), which this case bound. For the low 


value 11.1 the impact velocity the lower bound for 


has error and the upper bound has error 
The time which the moment becomes elastic again 

(t, 0.375) approximated 0.406 Eq. (59), derived 


from the lower bound, and 0.346 Eq. (56), derived 


from the upper seen that for low values 
the velocity the elasto-plastic theory needed pre- 
dict the behavior the beam. 
Fig. gives the elasto-plastic moment diagram between 
case the high value the initial velocity would create 
elastic bending moment times the value the capacity 
moment. 
Fig. gives the plastic angle together with its upper 


and lower bounds, which for the high value 67.4 the initial 
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= 


velocity parameter are identical within the accuracy the 
graphs except for the interval t=0, For such 


the value obtained from the bounds 


coincide with the true value, predicted the theory. 


values 


thus seen that the case impulsive motion the 
approximate solutions based upon rigid-plastic behavior become 
whenever the maximum theoretical bending moment due 


the impact much larger, say more times larger, than 


the capacity moment, while elasto-plastic solutions are re- 


quired all other cases. 


bs 
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ELASTIC 
STATES STRESS 


FREE-FREE BEAM DYNAMIC FORCES 
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